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ANALYSIS OF MICROWAVE 


RADIOMETRIC MEASUREMENTS FROM SKYLAB 

Robert M. Lerner and James P. Hollinger 
Naval Research Laboratory 
Washington, D. C. 20375 

INTRODUCTION: 

The earth-orbiting satellite SKYLAB was operated during 
the summer, fall and winter of 1973-74. Its orbit criss- 
crossed the earth at an altitude of approximately 235 nautical 
| i'"j ' : 

miles between latitudes SON and 50S. A group of five remote 

earth sensing instruments was designated the Earth Resources 
Experimental Package (EREP) . The results from two of these, 
S-193 and S-194, are the subject of this microwave sensor 
report. Their locations on SKYLAB are shown in Figure 1. 

The S-193 microwave sensor consisted of a scatterometer , 
altimeter and passive radiometer, all sharing a common 

parabolic dish antenna and operating at a frequency of 

■\! 

13.9 GHz. The antenna can be tilted 48° off-nadir in either 

- ■ j ' ' ■ ' ■ ' • ’ : . ' ' . . ■ ■ • • ' - . 

direction along the spacecraft ground track or across this 
track. The S-194 radiometer employs a fixed, one meter 
phased array antenna with its field of view centered on 
nadir. The operating frequency of the S-194 is 1.4 GHz. 

The objective of the investigation reported here is to 
establish the degree to which quantitative measurements of 
sea surface conditions and related wind fields can be made 

1 



using the S-193 and S-194 radiometers. The first priority 
qf the investigation is to demonstrate the reliable operation 
of the radiometers and the feasibility of satellite-based 

1 : ’• ' ‘ " ■ . i 

observations under ideal ocean and atmospheric conditions. 

The second priority is to examine the effects of marine wind 
speed, sea surface roughness and foam coverage on radiometer 
antenna temperature in order to evaluate the potential for 
passive microwave determination of ocean surface wind fields 
from a satellite. The third priority is to examine the 
effects of changes in sea salinity of selected target areas 
in order to test the capability of S-194 to determine sea 


salinity. The fourth priority is to examine the effects of 
a wide variety of cloud and precipitation conditions over 
the ocean in order tq investigate the effects of these 
atmospheric conditions on passive microwave measurements 
from a satellite. 

The microwave radiometer measures the absolute level of 
the total power received by its antenna over a particular 
frequency interval Af. Although the received power may 
contain spurious manmade radiations at times, we are concerned 
only with the absolute measurement of the self radiation of 
the various materials composing the earth’s environment. 

The signal received by an earth viewing radiometer is 
composed of not only the radiation emitted by the earth’s 
surface, but also the downwelling sky radiation reflected 
from the surface, as well as the emission and attenuation of 
the atmosphere between the surface and the radiometer. In 
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general, high gain antennas receive power primarily over a 
relatively narrow beam or solid angle, however some power is 
usually received from all directions. The power received in 
these side lobes and back lobes, though possibly small, must 

be taken into account if the intensity and spatial distribu- 

| ' • ' ‘ 

tion of the incident signal is to be accurately inferred 

from the measured total power. 

In general the measured intensity distribtion is the 
tirue intensity distribution convolved with the antenna 
reception pattern. There are two approaches to the inter- 
pretation of the radiometric measurements. One is to obtain 
the true brightness distribution from the measured distri- 
bution using a detailed knowledge of the antenna pattern and 
any one of a number of inversion techniques (e.g. 1). It 
should be noted that the restoration process can only provide 
ah intensity distribution which satisfies the convolution 
equation and the measurements but which is not unique since 
an infinite number of solutions can be found for any finite 
antenna pattern and set of measurements. Generally, the 
simplest solution consistent with the measurements is selected, 
or a priori knowledge of the intensity distribution is used 
to obtain an acceptable solution. An alternate approach, 
and the one adopted for this investigation, is to model the 
earth environment, calculate an intensity distribution, 
convolve it with the known antenna pattern, and compare the 
calculated received power with the actual measurements. The 
model is adjusted until the calculations agree with the 



measurements to within the measurement error. Both approaches 
are equally valid and general. The latter approach may be 
preferred if the measurements are made only at a few restricted 
positions. Since, in the case of S-194 measurements are 
made only at nadir, the inversion technique is not practical 
and the latter approach was selected for all of the studies 
reported here. 

The approach is to calculate the microwave intensity 
using theoretical models for the atmospheric radiation and 
attenuation, and for the ocean surface emission and reflection; 
and then to convolve this total radiation with the detailed 
antenna pattern. These calculations are then compared withj 
the measurements. In order to perform the convolution, it 
is necessary to keep track of the polarization and intensity 
of the various radiations from the environment and properly 
reference them to the polarization and gain of the antenna. 

This is conveniently done by using Stokes parameters (e.g. 

2) and translating from the earth oriented coordinate 
system, with respect to which the environmental calculations 
are made, to a satellite-borne antenna coordinate system. 

These calculations will be described in detail and then 
applied to measurements of increasingly complex environmental 
conditions as follows: 

1. Use the model and antenna pattern to correctly 
calculate recorded antenna temperature over ideal ocean 
conditions, temperature over ideal ocean conditions. 

2. Make beam corrections in the vicinity of rapidly 
changing signal such as land-Water and water-land interfaces. 
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3 . Measure the effect of changes in the salinity of 
the target area. 

4 . Examine the effects of changing thermodynamic water 
temperature. 

5. Examine the effects of sea surface roughness and 
foam coverage. 

The analysis of this report will be confined to measure- 
ments from S-194. Unfortunately, the detailed, accurate 
antenna pattern necessary for the analysis of S-193 data was 
not received until the end of this investigation and after 

i 

funding had been expended. However, the complete environmental 
model and formalism necessary for the analysis of S-193 data 
hs been completed and is included in this report. It is 
hoped that analysis of the S-193 data can be accomplished at 
some future date using alternate funding. 

SENSOR DESCRIPTION: 

The S-193, 13.9 GHz radiometer shared a 48-inch diameter 
parabolic antenna with an active microwave altimeter and 
scjatterometer . It was capable of receiving both horizontally 
and vertically polarized radiation. The half -power beamwidth 
of the S-193 was approximately 1.6°, equivalent to a nadir- 
locking, half-power footprint of approximately 6 nautical 
miles assuming a satellite altitude of 235 n.mi. The antenna 
is mounted in such a fashion as to make possible off-nadir 
scans up to 48° corresponding to 50.3® with respect to local 
vertical at the earth's surface. 
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The available scan modes were in-track contiguous 
(ITC) , cross-tract contiguous (CTC) , in-track non-contiguous 
j(ITNC) , cross-track non-contiguous right (CTNC-R) , cross- 
track non-contiguous left (CTNC-L) , and cross-track non- 
contiguous left and right (CTNC-L/R) j The S-193 Sensor 
Performance Evaluation Group (S-193 SPEG) (3) determined 

thkt the standard deviation for measurements using the S-193 

■ 1 . .. : ■ 

(radiometric precision) had an average value of 1.62°K. 

'. | ' • '' ; . . 

Dynamic range of the radiometer receiver was determined by 

j i : . ' I • 

the S-193 SPEG as being 4 to 353°K with a linearity within 

j ’ . ■ : ] . _• • ■ . ■ 

4.3°K of best straight line over 350°K. The footprints for 
S-193 in the ITNC and CTNC-R modes are shown in Figure 2. 

The S-194 operates at a center frequency of 1.414 GHz 

• I 

with a bandwidth of 27 MHz. It utilizes a square, planar, 
8x8 element phased array antenna measuring one meter on a 
side. The half-power beamwidth was 15° with greater than 
97% main-beam efficiency. This results in a half-power 

• • ’ ' j I ■ ; -T \ \ : 

footprint of approximately 60 n.mi. The S-194 footprint is j 

• i ; ! 

I . ■ ; : _ i 

represented by the large circle in Figure 2. A. I.L. provided 

i ; ■ 

i . 

us with very complete antenna patterns taken before launch 

of the flight version of the S-194 antenna. The patterns 

■ ' . j"''' . ■' • ■ . .. 

consist of 72 cuts through the main-beam axis spaced every 5 

• ' : • j • 

degrees in azimuth. The normalized S-194 antenna response 
pattern, integrated over 360° in azimuth, as a function of 
the angle off the main-beam axis, theta, is presented in 
Figure 3. Also given is the percentage effective solid 
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angle contained within theta as a function of theta. Although 

j ! " ' • •' j ’ ! ’ : 

the pattern averaged in azimuth is given in Figure 3, all of 
the detailed pattern cuts were used in the analysis of S-194 
data. The radiometer utilizes a calibration scheme referenced 
to a fixed "hot" and "cold" load input. The data are pulse- 
code modulated and recorded on magnetic tape. The "hot" and 
"cold" reference noise generators provide inflight calibration 
of the radiometer electronics. The accuracy and precision of 
the radiometric measurements depend on the stability and 
accuracy of these calibration sources. 

S-194 CALIBRATION : 


The first deep-space calibration pass during SL-2 
provided a measured antenna temperature greatly in excess of 
the known L-band cosmic background temperature. The S-194 
Sensor Performance Evaluation Group (S-194 SPEG) (4) determined 
that the pre-flight electronics box calibration must be in 
error. As a matter of paracticality / it was assumed that 
all insertion loss terms remained unchanged from their 
prelaunch measured values, and that the total error in the 
deep-space calibration could be assigned to the electronics 
box calibration constants C 2 and C^. This was done since 
ail the loss terms including C 2 and are interdependent 

■ . j ■ i ; 

and no one term can be determined uniquely from the flight 
data available. The system processing algorithms were 
reworked using a known value for cosmic background temperature 
and some pre-flight hot calibration source temperatures. 
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Figure 3 


PERCENTAGE SOLID ANGLE WITHIN THETA 



Simultaneous solutions of the reworked system equations 
yielded new values for C 2 and C^. These new values were 
applied to the data processing for all three Skylab missions 
before the data was released. 

As considerable amounts of data from SL-2, SL-3, and 
SL-4 became available, it was obvious that a second correction 
iteration was necessary. C 2 and C 4 were again adjusted 
using a slightly better value for the cosmic background 
temperature and a value from Pass 56 for the antenna tempera- 
ture of the Gulf of Mexico. This adjustment was not included 
in the released data. Thus it was necessary to change each 
output value individually to incorporate the new values of 
C 2 and in the antenna temperature determination. 

After reducing the available data using these values 
for C 2 and C^, provided by the S-194 SPEC, and comparing 
with our calculations we noted a systematic offset of 
approximately 3.6 degrees. We then elected to further 
refine the system equations by using the values of 2.8°K for 
cosmic background temperature (18) and 95°K for the Gulf of 
Mexico temperature. This calculation resulted in what we 
now consider the best available values for and . 

Utilizing these values provides for close agreement between 
all measured and calculated values. It may be argued that 
the adjustment in the calibration constants based upon the 
calculated value for the temperature of Gulf of Mexico is 
somewhat of a bootstrap operation and subsequent agreement 
between the recalibrated measurements and the calculations 
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is therefore not surprising. We feel this procedure is 
appropriate and analysis of the data is valid since the five 

j ^ 

measurements over the Gulf of Mexico used in the recalibra- 

j 

\ 

tion constitute a very small set of the total measurements 
analyzed. Further, they were taken under nearly ideal 
conditions of clear weather and calm seas and excellent in 
situ ground truth is available from ships and aircraft. 

i 

Values of C 2 and were computed for each of the five Gulf 
of Mexico passes occurring on 6/10, 6/11, 8/8, 1/6, 1/8 
1973. The average values of C 2 = 6.2124, = 13.0694 which 

were deduced were then used to reduce all the available 
data. Due to the significance of these changes in the 
calibration constants, we feel it is worthwhile to summarize 
the description of the procedure here even though it appears 
in S-194 SPEG report (4). 

The various temperature and .losses necessary for the 
calibration analysis are shown in Figure 4 which is similar 
to Figure 3 . 2 .1.3-1 of MSC-05528. The hot reference noise 

■''j •: ' ; ' ' 

generator (RNG) was an accurately known and highly stable 

J 

temperature controlled device designed to provide a reliable 

■ j • • " 

i \ 

reference for evaluating the stability of the radiometer 

system. In order to check the accuracy of the hot RNG, a 
comparison was made before Skylab flight between the true 
radiometric temperature of the source, and the Y-f actor 
output from the radiometer obtained during the hot calibra- 
tion mode, corrected for cable losses. The relationship 
between the radiometric temperature at the diode switch and 
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the computed source temperature is given by the following 
expressions. 


T = 
SH 


t ph ~ ^ 1 ~ L CH^ 
l ch 


( 1 ) 


T PA T SA (L GA ) + (1 " L CA ) fc CA 


( 2 ) 


m — «T» 

a PH A PA 


(3) 


SH 


T SA (L CA* + (1 ~ L CA ) fc CA ~ ^ 1 ** L CH^ 

l ch 


(4) 


where 


^CH ^ T SH + fc D 


(5) 


1-L 


t = t + 
-GA L CAA 


CAB 


1-L 


(t D " t CAA ) 


( 6 ) 


CA 


Several basic assumptions were made in attempting to correlate 
the computed source temperature with the hot RNG source 
thermal temperature: 

1. The thermal temperature of the hot RNG is its true 
radiometric temperature. 

2. The loss factor for the cable which connects the 
hot RNG to the diode switch is well known. 

3. The temperature of the external RF cable ( T CAA ) is 
the same as the temperature of the S-194 antenna 
'structure. 
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4 . 


In the hot calibration mode T p ^ = T PA* T he s "^^ 4 


is calibrated such that for any radiometric tempera- 
ture applied to the input of the RF antenna calbe 

would be identically recorded in the output data. 

| 1 

A discrepancy of 1.35°K between the computed source 

temperature of 37l.49°K and the hot RGN thermal temperature 
of 372.84°K was recorded. This discrepancy was apparently 
due to inaccuracies in Equations (1) - (6) above. 

i • . ■ ; I ■ 

A study of the raw data for the hot calibration periods ! 
during SL-2, SL-3, and SL-4 gave a value for the standard 
deviation of the radiometer output TSA of a ^ = 0.1329°K, 
a 3 = 0.1538°K, and a ^ = 0.1947°K over all the calibration 
passes for which the hot source has stabilized. There is 
currently no explanation for the slight increase of a from 
SL-2 to SL-4. 

I A calibration curve relating the radiometric temperatures 
present at the Dicke switch and the corresponding output Y- 
f actors can be computed from the equation* 


a a a c 

T -T 
A PA A PC 


^ - A c 

T -T 
PH PC 


(7) 


Solving for T pA 


A A ~ A C 

T = T + — — (T 

A PA PC An - A c U PH 


T PC^ 


(8) 


*See reference Table I for definition of terms, (p.26) 
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The relation between the hot RNG radiometric temperature 


at the Dicke switch (T_.„) and the hot calibration radiometric 

rn 

temperature (T must include cable losses. 

Oil 


t ph t sh l ch + (1 l ch ) 


(9) 


This is equally true for the cold calibration radiometric 
temperature T gc , 


T PC T SC L CC + ^ L CC^ fc CC + “ T SC^ ^ “ L CCB^ ‘ 


Similarly, 


T PA T SA L CA + (1 “ L CA ) t Ck * 


(ID 


There exists a relationship between the - RF cable thermal 
temperature t CAA and the orbital Beta angle ( 3 ) which is the 

I • 

I . ' . 

angle between the sun vector and its projection into the 
ofbital plane. This relationship was determined emperically 
for a t Q = 297°K and for Beta angles between zero and 53 


degrees. 


'CAA 


- [** P - 10 ][o.86 - %§i^ + 50 , 


where 


tg = thermal temperature values in word 9 of data frame , 
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tiQ = thermal temperature value in word 10 of data frame 


Substituting for T pA , T pc , and T pH in equation (8) using 
equations (9) , (10) , and (11) yields 


SA CA + (1-L ca ) t CA T sc L CC + (1-L CC ) t cc + (t 0 -Tg C ) 


a a' a c 


(1 “ L CCB ) + A h “A c [ t sh l ch + ( 1 “ L CH ) fc CH 


_T SC L CC “ {1 “ L CC ) ” ^ t o _T SC^ ^” l ccb^J 


for t sa / 


SA 


” L ca [ T SC L CC + ^^CC^CC ” ^ 1-L CA^ fc CA + ^ t D" T SC ) {1 “ L CCB^. 


+ L CA (A -a ) [ T SH L CH T SC L CC + ^ 1 ” L CH^ fc CH 


(l-Lcc)tcc “ ^“^CCB^ ] * 


Now define: 


r - L cc 
1 ' ^ ’ 


C 2 5 L 


t [ (1 - L cc )t cc 


(1 “ L CA* fc CA + (t D" T SC ) * 1_L CCB*_ / 


c, = 


, CH 


L CA ’ 
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(12) 


C 4 E 


= r u- 

l ca L 


l ch )1: ch (1 " l cc )1: cc ' (t D~ T SC ) (1 ” l ccb ) 


Hence, (Eq. 12) may now be written as 


t sa c i t sc + c 2 + a h -a c " c i T sc + c a \ 


( 13 ) 


Correction terms reflecting the losses in the RF cable and 
the changes in radiometer performance resulting from changes 
in operating temperature must be added to equation (13) . 


? SA C 1 T SC + C 2 + A h “A c (C 3 T SH ' C 1 T SC + C 4 } 


( 1_ L CA } ( t C AC AL” t CA * 


(14) 


CA 


where 

t CA cAL =t W at the t: ^ me dee P space data acquisition 
- 295°K. 

Utilizing the following values from EREP 11 lunar calibration 
pass, it is possible to solve for T gA : 


T SH = 373 . 0875°K, A fi = 15 , 
T gc = 249 . 0179°K, A c = 170, 
t D = 288 . 4948°K, I 

t^ = 233. 5538°K, 

^10 = 236. 9276°K/ 

T sa * 382.9161 - 0.8207 A A . 
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At orbital midnight, A A = 420. Therefore at orbital 
midnight T gA = 38.24°K. 

The antenna corrected temperature (T AC ) is as follows: 


T SA ~ (1 " L A) t L 


AC = 


jvif ; 

L r a +1 J l 


Vt* + - (1 


-v**] 


i .felil 

l r a + i J J 


where 

i\ 

= 6 (1,5t W9 + 1 * 5t W10 + ^5 + *~W6 + fc w8 * 

See Reference Table 1 for definitions of terms. 
Therefore, in terms of T gA , 


‘SA 


“AC 


0. 911’ 


- 23.40, 


for T sa = 38 . 24°K, T AC + 18.54°K. 

This is not an acceptable temperature for the cosmic 
background radiation which is known to be approximately 
2. 8°K at 1.4 GHz (18). S-194 SPEC felt that the largest 

part of this error was due to inaccuracies in the electronics 
box calibration. Re-calibration was first done using 4°K as 
deep space brighness temperature and data from the AIL pre- 
installation calibration tests for the hot source . 
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(15) 



t sa C 1 T SC + 


C 2 + A c -A r [ C 3 T SH " C 1 T SC + C 4] 


A C~ A H [ C3TsH ~ C 1 T SC + c 4]‘ 


Let 


T k C 1 T SC + C 2 ' • 


K 


a c” a h t 


C 3 T SH " C 1 T SC + C 4 * 


It follows that 


t sa - (t k + A c K > - A A K • 


*Now taking the hot and cold reference temperatures and the 
corresponding Y- factor values for these temperatures T R/ K, 
Cj, and C 4 may be solved for. 

m 

SAH is the external hot reference temperature and A AH 
is the corresponding Y-factor value. 

T SAC the externa l c °l d reference temperture and a ac 
is the corresponding Y-factor value. 

’■ ; ■ 1 ... . I, .. 

T SAH = (T k + A C K) - A ah K 


SAC 


’k + A c K ) a ac k 


Prom the pre-installation test calibration 

T„ah = 372. 2°K and A AH = 13.2 bit counts. 


( 16 ) 

(17) 

(18) 


(19) 

(20) 
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Assuming a deep space temperature of 4°K 


T AC 0.9117 " 23 . 40°K - 4°K . 

T SAC = (23.40°K + 4°K) 0.9117 = 24.9806°K 
A ac =420 bit counts. 


Solving simultaneously (19) and (20) gives 
T R = 238.3547°K; K = 0.8536°K per bit. 

.1 ■ ’ ■ ' ■ ; i 

K is a measure of the sensitivity of the radiometer. 

Applying these values to equations (17) and (18) and using 
the values 

T SH = 273 . 0875°K, A R = 15 

T ge = 249.0179°K, A c = 170, = 0.9615 and C 3 = 0.9725; 

and C 4 are obtained. 

C 2 = -1.0749, 

C 4 = 8. 9017. :■ 


As a result of recalibration, the correction for C 2 due to 

antenna cable loss temperature becomes (270.2 - t^) 

These corrected values for and C. derived from an assumed 

2 4 

deep space radiation temperature of approximately 4 °K hopefully 
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were to result in measured antenna temperatures close to 
those expected from modeling and previous aircraft observa- 
tion. Unfortunately examination of data from SL-2, SL-3, 
and SL-4 which were processed with these values of C 2 and 
C 4 showed discrepancies of approximately 10° K from expected 
values of antenna temperature. 

A second adjustment to the "constants" C 2 and was 
then computed using data from Pass 56 over the Gulf of 
Mexico. This second recalibration was performed using the 
same general method as was used for the first recalibration 
try after SL-2. The major difference was that this second 
recalibration used only flight data rather than the combina- 
tion of flight and test data used for the previous recalibration. 

i 

Using the transfer equation (15) for antenna temperature 


(T ac ) / temperatures at the radiometer antenna input cable 

i • _ • 

(T sa ) were computed from two calibration targets: Deep 

Space Pass LC-5 (T AC = 3.3°K) and the Gulf of Mexico Pass 56 
(T ac = 99. 0°K) . Upon solving the equations for T gA , the 
results were: 


T SA = 23.90°K for Deep Space Pass LC-5 


T SA = 111. 12°K for Gulf of Mexico Pass 56 


Using these two values of T gA and the equation 


T SA = C 1 T SC + C 2 + (C 3 T SH - C 1 T SC + C 4> 

(*c - a h> V : 
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17.995. 


yields solutions for C 2 and C^, C 2 = 14.425, = 

A new correction factor to C 2 from the cable loss 
temperature t CA must also be used 




(252.5 - t CA ) 



To verify the accuracy of these two re-recomputed 

•i . . : } 

calibration constants, they were applied to six selected 
targets taken from SL-3. For these six targets using the 
existing calibration constants (i.e. C 2 and C 4 after the 
first recomputation which was incorporated in TR.«524 outputs) 
the discrepancy between the measured antenna temperature and 
the S-194 SPEG modeled computed antenna temperature averaged 
approximately six percent. Using the new values for C 2 and 
C 4 the discrepancy is reduced to approximately seven tenths 
of a percent. A least square fit line for the SL-3 selected 

j . • ' 

sites shows good linearity over the total dynamic range when 
the least values of C 2 and C 4 are used. 

The computer-compatible tapes, tabulations and plots 
sent out to the Principal Investigators for S-194 are in 
error in that they are not computed using the latest values 
for C 2 and C 4 calculated by the S-194 SPEG . Each investi- 
gator must thus change each output value individually. 

Upon utilizing these latest calibration constants in 
the reduction of the available data, . a consistent offset of 
approximately 3.6°K was noted between the "corrected" data 
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and our calculations. Analysis of the method used in the 
recalibrations led to the substitution of a lower and more 


realistic radiometric antenna temperature for the Gulf of 
Mexico. Details of the target model used by the S-194 SPEG 
are not presented in MSC-05528. Our model, using data from 
Ho et al. (6) that may not have been available to the S-194 
SPEG, provides an antenna temperature of 95°K for the Gulf 
of Mexico. This value is four degrees less than the 99°K 
used by the S-194 SPEG. 

In order to calculate a correction algorithm for the 


corrected antenna temperature (I^J , it is necessary to 
combine equations (14) and (15). Correct T gA , T AC , and 
C 4 are defined as T gA , Tj^ c , C^, . 



Substituting the following values: 

C x = 0.9615 
CJ = $*21.24 
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C 2 = -1*0749 


C 3 = 0.9725 

= 13.0694 

C 4 = 8.9017 

L- = 0.9231 
A 

R. = 1.25 
A 

yields 





TABLE 1 


SYMBOL 


DEFINITION 


J CA 


CAB 


'CA 


'CAA 


'CH 


Sj5' S?6' fc W8 ' 

w9 ' ^WIO 


CH 


‘PA 


PH 


Electronics Box Calibration Constants, 
Evaluated from Cable Loss Data 

Electronics Box Calibration Constants, 
Evaluated from System Measurements 

Loss Factor of Cables Connecting Antenna 
Terminal to Diode Switch 

Loss Factor of Portion of Antenna RF Cable 
Within Electronics Box 

j • . 4 ' . 

Loss Factor of Antenna 

j . 

VSWR of Antenna at its Output Terminal 

Average Temperature of Cable Connecting 
Antenna Terminal to Diode Switch 

Temperature of External Portion of Cable 
Connecting to Antenna Terminal 

Electronics Box Enclosure Temperature 

Average Temperature of Cable Connecting Hot 
RNG Source to Diode Switch 

■ • .j ' ' ' 

Effective Loss Temperature of the Antenna i 
Five Antenna Thermal Temperature Monitors 


Loss Factor of Cable Connecting Hot RNG 
Source to Diode Switch 

Antenna Terminal Radiometric Temperature 
Referenced to the Diode Switch 

Hot RNG Radiometric Temperature Referenced 
to Diode Switch 


SH 


PC 


SC 


Hot RNG Source Temperature 

Cold RNG Radiometric Temperature Referenced 
to the Diode Switch 

Cold RNG Source Temperature 


VALUE 


0.9615 

0.9725 

“ 1 . 0749 * 

8 . 9017 * 

0.97 

0.99 

0.9231 

1.25 
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TABLE 1 (continued) 



SYMBOL DEFINITION 


VALUE 


T 

SA 

T AC 

a A 

A C 

L cc 

l ccb 


Antenna Terminal Radiometric Temperature 
Referenced to the Input of the RF Antenna 
Cable 

Corrected Antenna Radiometric Temperature 

j 

Antenna Radiometric Measurement 
{Y-Factor Bit Count in Mode 1) 

Hot Calibration Radiometric Measurement 
(Y-Factor Bit Count in Mode 2) 

Cold Calibration Radiometric Measurement 
(Y-Factor Bit Count in Mode 3) 

Loss Factor of Cable from Cold Source to 
Cold Source Box Terminal 

Loss Factor of Cable from Cold Source Box 
Terminal to Diode Switch 

Thermal Temperature of Cable from Cold 
Source to Cold Source Box Terminal 

* These values are used in the TR-524 

outputs. The values which are now 

considered to be more correct are 

C, = €.2124 

CZ = 13.0694 
4 




DESCRIPTION OF CALCULATIONS: 

A computer program ANTTEMP was written to perform the 

| " 

calculations necessary to model the total brightness distribu- 
tion seen by the antenna and convolve this brightness distribu- 
tion with the antenna pattern in order to evaluate Skylab 
antenna temperature measurements. The calculations necessary 
to model the brightness temperature are as follows: (1) the 

radiation emitted by the sea, (2) the upwelling radiation 
emitted by the atmosphere, (3) the downwelling radiation 
emitted by the atmosphere and extraterrestrial sources and 
reflected upward from the sea surface, (4) the attenuation 
of the radiation emitted and reflected from the sea surface 
. by the intervening atmosphere, (5) the transformation of the 

'. 'I ■ • . ■ ..... ; • ' . i '• 

antenna properties and orientation from the satellite 
position and reference frame to the coordinate frame on the 
sea surface with respect to which (1), (2), (3), and (4) are 

•'■I' ■ . ■ 

most easily calculated, and (6) the convolution of the 
antenna reception pattern with the total radiation at each 
position over the earth's surface. Each of these calcula- 
tions will be described in turn. 

In order to calculate the emission and reflection of 
microwave radiation by the sea surface, the dielectric 
constant of sea water must be known. This is calculated by 
a subroutine DIELCON which requires 'the inputs of sea 
temperature and salinity and the observational frequency. 

The complex dielectric constant of water e may be written 
as> .. 
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E a £ o (e* - je") 


where £ is the dielectric constant of a vacuum. The real 
o 

and imaginary parts of the relative dielectric constant of 
water i s' and e” , are well represented by 

, < E S - 

c . = —2 + 

1 + (27TfT)^ 

and 


(e_ - £ ) (2-irfi) 2a 

— * + — (22) 

1 + (2TTfl)^ f 


where f is the frequency in hertz, t is the relaxation time 
.■ j ; - ‘ j 

in seconds, a is the ionic conductivity in ESU (1 ESU = 4 xe q 

MHOS/meter in MKS units) , e is the static relative dielectric 

• j ■ ■ ■ ' ' . 

constant and £ is the relative dielectric constant at 

| CO 

very high frequencies. Both £ and £ are dimensionless. 

S 00 

The constant value of 4.9 given by Saxton and Lane (5) 
is used for £ . 

CO 

The static dielectric constant, relaxation time, and 

... ....... 

conductivity are functions of salinity and temperature. To 
calculate the dielectric properties of water for any combi- 

.i • ' ■ j 

riation of salinity and temperature, a second order polynomial 
ijn salinity and temperature was fitted by the method of 
least squares to the experimentally determined values of e , 

T t and a given by Saxton and Lane (5) . 
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This determination of the dielectric properties of 
water as a function of salinity and temperature is probably 

s 

good to better than one percent at the frequencies of S-193 
and S-194. However even greater accuracy is possible at 
1.43 GHz tanks to the recent accurate measurements of Ho, 
Love, and Van Melle(6). Their measured values are accurate 
to better than 0.2 percent in s' and 0.4 percent in e". 

Therefore the above determination is used for analysis of 

| 

S-193 data but for the analysis of S-194 data we use the 
results of Ho et al. They give the dielectric properties of 
sea water at' 1.43 GHz as 

e' (x,T) = (e^ + a Q +.a 1 x-l)/(a 0 + a^) 

e" (x,T = (C Q + C X X) Ujj - l)/(a Q + a^x). (23) 

where s' is the real part of the dielectric constant for 
w 

distilled water at sea temperature T in degrees Celsuis; X 
is the chlorinity in PPP which is related to the salinity S 
by X - (S - 0. 03)/l. 805; a Q , a^, Cq, and are coefficients; 
and 

£ ' 

e w 

£*• 

W 
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= 85.98 - 0.271 T - 3.70 x 10“ 3 T 2 + 6.0 x 10 5 T 3 
= 13.03 - 0.5062 T + 9. 91 x 10“ 3 T 2 - 7.3 x 10' 5 T 3 



a Q = 1.0022 

a x = 0.005786 - 1.96 x 10" 5 T 

C Q = 0.1564 - 4.12 x 10" 3 T + 2.07 x 10" 5 T 2 + 5.13 x 10 _7 T 3 
C x = 0.02231 + 1.105 x 10“ 3 T - 9.63 x 10‘ 6 T 2 + 4.18 x 10“ 7 T 3 . 


Using the results of subroutine DIELCON and the angle 
of incidence on the sea surface 0, the Fresnel reflection 
coefficients for a specular air-water interface are given by 



(cos 9 - P) 2 + Q 2 
(cos 0 + P) 2 + Q 2 




(e 1 cos 9 - P) 2 4- (e" cos 6 - Q) 2 
(e ' cos 9 + P) 2 + (£" cos 0 + Q) 2 


(24) 


The subscripts h and v refer to the horizontal and vertical 

lijnearly polarized components of the radiation defined as 

. j ' • . ■ 

bejing when the electric field vector is perpendicular or 

i ; 

parallel to the plane of incidence. The variables P and Q 
are given by: 


P - sin 2 0) 2 + e" 2 ' + (e 


and 

Q = 1 1 [^(£' - sin 2 6) 2 + e " 2 - ( £ 1 


sin 2 0)J | 


sin 2 0) J | 


1/2 


1/2 
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Tabulated and plotted values of the dielectric constant and 
the Fresnel reflection coefficients, using the dielectric 
properties based upon the measurements of Saxton and Lane, 
have been calculated from the above equations as a function 
of temperature, salinity, look angle and frequency and are 
available in a separate' report (7) . Calculations of the 
brightness temperature of a specular sea, as a function of 
sea salinity and sea temperature at 1.4 GHz, based upon the 
dielectric properties measured by Ho, et al., are given in 
Figures 5 and 6. These calculations do not include atmo- 
spheric effects and are intended to display the quantitative 
dependence of the brightness temperature on sea salinity and, 
sea temperature for S-194. Note in particular the near 
independence of the brightness temperature on sea temperature 
for salinities near 35 parts per thousand; typical values of 
ocean salinity. 

In order to sum up the various radiations and convolve 
the sum with the antenna pattern, it is necessary to keep 
track of the polarization and intensity of the radiations. 
This is conveniently done by specifying the radiation in 
terms of Stokes parameters since the Stokes parameters for 
the sum of several independent radiations is the sum of the 
Stokes parameters of the separate radiations. The Stokes 
parameters of the radiation emitted from a specular sea are 
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BRIGHTNESS TEMPERATURE ( £ 'K) 



Figure 5 




SEA TEMPERATURE (°C) 


Figure 6 
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( 25 ) 


I bb = (2hf/X 2 )/(e hf/kT - 1) . (26) 

Here h and k are Plank's and Boltzmann's constants respec- 
tively, X is the wavelength and T the sea temperature. In 
most cases of microwave radiometry of the earth, and parti- 
cularly at the frequencies of S-193 and S-194 the Rayleigh- 
Jeans approximation to Plank's blackbody radiation law is 
sufficiently accurate, especially when, as described later, 
the appropriate atmospheric boundary condition is used. We 
have then, 

I hh = 2kT/X 2 . (27) 

Note that the radiation from any source can always be 
described in terms of an equivalent blackbody temperature or 
brightness temperature. The brightness temperature T g is 
the temperature to which a blackbody must be raised in order 
that the blackbody radiates the Same intensity, at frequency 

; -3a 




i 

£, as does the source. This is particularly convenient when 

I ...... 

the Rayleigh- Jeans approximation is valid, for then the 

i 

intensity is directly proportional to the brightness tempera- 
ture . 


The Stokes parameters of the downwelling atmospheric 


radiation, Ig^y* reflected from the sea surface are given by 


II 

n 

H 

kr£ + r;)/2] 

■^sky 

°r = 

ItR h - R v ,/2) 

I sky 

u = 

0 


r 


ii 

u 

> 

0 



The Stokes parameters for the emission and reflection 
ojf radiation by the sea surface given by the above equations 
refer to a specular sea. The effects of ocean roughness 
require a modification. 

The microwave brightness temperature dependence on 
surface roughness (wind speed) arises from two effects. The 
first effect results from the increasing roughness of the 
compact water surface and the second effect from the increasing 
coverage of white caps and sea foam with increasing wind 
speed. Both effects increase with wind speed and frequency 
and depend upon the incidence angle and polarization. The 
Surface roughness effect is closely coupled to the local 
wind field; rapidly responding to changes in the local wind 
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and hence is relatively insensitive to the energy content of 
low frequency gravity waves. Evidence indicates that it is 
primarily dependent upon the mean square surface slope and 
is thus due to surface waves of dimensions comparable to and 
shorter than the observational wavelength. The sea foam 
effect results from the relatively high microwave brightness 
temperature of sea foam and its increasing coverage with 

wind speed. The determination of the brightness temperature 

! 

dependence on wind speed is complicated by the fact that the 
foam coverage of the sea surface and the roughness effect 
depend, not only on the local wind, but also on the air-sea 
temperature difference, the duration and fetch of the wind, 
as well as on the history of the wave spectrum of the sea 
area being observed. 

No completely satisfactory analytical model for the 
microwave brightness temperature dependence on ocean rough- 
ness is available. Therefore the empirical expression 

AT B = 0.134 AW f 1/2 (29) 

was derived from measurements. It relates the change in 
brightness temperature to wind speed as a function of 
frequency. Here AT_ is the increase in brightness tempera- 
ture in °K, AW is the change in wind speed in knots and f is 
the observational frequency in GHz. The measurements upon 
which this expression is based include observations from a 
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bridge, an ocean tower and seyeral different aircraft. They 
were made at varying incidence angles and polarizations , 
over a wide range in atmospheric conditions and using a 
Variety of antenna and calibration procedures and show 
considerable scatter. Thus the above expression is hardly 
exact. But it does provide a reasonable estimate of the 
roughness effect at 1.41 GHz. The wind speed, dependence 

derived from it at 1.41 GHz of 0.16°K/knot is in good 

• | ’ . ! 

agreement with measurements from an ocean tower (8) and a 

bridge over Cape Cod Canal (9) . Since the wind speed effect 
changes little with incidence angle out to angles off nadir 
well in excess of those for which the S-194 antenna can 

I l : i 

receive significant signal, no incidence angle dependence 
was incorporated. The wind speed effect was taken into 
account in the analysis of the S-194 data by subtracting the 
ratio of AT divided by the sea temperature from the Fresnel 
reflection coefficients calculated for a specular sea. 

As a first attempt to include the wind speed effect in 
the analysis of S-193 data, the geometric optics model as 
developed by Stogryn (10) was adopted. A subroutine STOGRYN 
was written to calculate the Stokes parameters for emission 
and reflection. Since no S-193 data has been analyzed, it 
has not been used. 

The downwelling and upwelling radiation emitted by the 
atmosphere as well as the attenuation of radiation by the 
atmosphere are calculated in a subroutine ATMOS. There are 
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three basic parts to ATMOS, one in which the calculation of 
absorption properties of significantly absorbing atmospheric 
constituents as a function of density, frequency, temperature, 
and pressure is done, a second part which distributes these 
parameters according to ah atmospheric model , and a third 
which relates the resulting distribution of physical tempera- 
ture and absorptivity to a radiometric temperature at any 
microwave frequency. 

For an atmosphere in local thermodynamic equilibrium 
and wherein negligible microwave scattering occurs, the 
radiometric properties may be described by an absorption 
coefficient which gives the amount of energy absorbed per 
unit of length along a particular direction in nepers per 
kilometer. This absorption coefficient is primarily dependent 
on oxygen, water vapor, and liquid water in the atmosphere 
and has been expressed as a function of the densities of 
these constituents and the ambient temperature and pressure 
for use in ATMOS. 

The absorption formulas are computed in three separate 
subroutines: oxygen, water vapor, and liquid water. 

Oxygen 

The calculation of microwave oxygen absorption is based 
on Meeks and Lilley's (11) extension of Van Vleck's (12) 
work with the modification suggested by Reber (13) . The 
model used in ATMOS takes into account 46 fine structure 
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transitions of the O 2 molecule, 45 of them from 48 to 71 GHz 
and one at 118 GHz. Pressure line broadening and temperature 
dependence are also accounted for. The basic formula relating 
the absorption coefficient a in nepers per kilometer to 
frequency v in gigahertz, pressure P in millimeters of 
mercury, and temperature T in degrees Kelvin is: 


CPv' 


a - 


2X «*p (—) • 

N w \ kT/ 


(30) 


Here N is the rotational quantum number which assumes only 
odd values from 1 to 45, 


"m n(n+1} 

= 2.06844 — ■■ m ■■ is the exponent in the Boltzman 
kT T 


distribution factor, 


S.. = F XT ,a M , + F. t a„ + F a XT is the line strength factor 
N N+ N+ N- N- o No 


Av 


Av 


N+ (V N+ " V)2 + Av2 (V N+ + V)2 + Av2 


F = 
* N- 


Av 


Av 


(v N _ -v) 2 + Av 2 (v N _ + V) 2 + Av 2 


V = -5 


Av 


v + Av 


2 ' 


v^ + is the frequency in gigahertz of a transition from 

total angular momentum quantum state (including spin 
of the two unpaired electrons) J = N to J = N + 1 , 
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v^_ is the frequency in gigahertz of a transition from 
total angular momentum quantum state J = N to 
J = N - 1, 

a = N(2N + 3) 

N+ N + 1 


_ (N + 1) { 2N - 1) 

a N- N 


2 (N + N + 1) (2N + 1) 

N(N + 1) 



Av = j 52.7 + 0.627 (P - 20.7) | x 10“ 3 GHz or 

c = 0.61576. 

Height, temperature, pressure, and frequency are inputs 
to the oxygen routine with the absorption coefficient as 
output. This analysis does not take into account high 
altitude Zeeman splitting (above 40 km) nor the doppler 
broadening above 80 km. However, since comparatively little 

•j ' 

absorption takes place due to low pressure at high altitudes, 

• " ■. j ‘ 

Zeeman and doppler may be neglected in almost all remote 
sensing applications except possibly high altitude upward 
looking radiometers for which a more involved computation is 
necessary. 


88 P 


too " 

T 


x 10 GHz whichever is less, and 


Water Vapor 


Microwave absorption due to water vapor is based on 
Staelin's (14) interpretation of Barrett and Chung's (15) 
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formula which in turn uses experimental results of several 
authors to extend Van Vlecks's (16) work on the subject. 

There is one principal rotational transition for water vapor 
which signigicantly affects microwaves and it has a frequency 
of 22.235 GHz. The equation used in ATMOS to calculate a in 
nepers per kilometer is: 





Liquid Water 

The absorption coefficient for liquid water in clouds 
is based on Paris's (17) Thesis. First the real and imaginary 
parts of the dielectric constant of water at a particular 
temperature and frequency are computed in a separate sub- 
routine which has a lower bound of -10®C, and then the 
absorption coefficient is computed as: 

P L v. (32) 

a is in nepers per kilometer, 

and E£ are the real and imaginary parts respectively 

of the dielectric constant of cloud water, 

■ i- ■ 3 

p T is the density of liquid water in gm/m , and 

Li 

v is the observation frequency. 

This approximate solution is valid for the microwave region, 
and non-raining clouds wherein the drop size is much smaller 
than the wavelength. (Rayleigh’s limiting case of dielectric 
spheres in a plane wave field.) The scattered component, it 
will be noted, is negligible and what remains is the 
component actually absorbed by the liquid water droplets. 
Temperature, frequency, and liquid water density are inputs 
to this subroutine and the absorption coefficient is the 
output . ■; 


a = 0.0629 


3E. 


(2 + 


V 


+ e: 
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The model atmosphere used in ATMOS has a linear tempera- 
ture lapse rate from the specified ground temperature to the 

specified tropopause height and temperature. The tropopause 

. j 

is : taken to be at a constant temperature. An exponential 
pressure decay derived directly from the hydrostatic equation 

. . i '■:■■■■■• ' 

is set at the ground by the local barometric pressure i 
Oxygen is assumed to be in constant proportion and thereby 
specified by the pressure vs. height array. Water vapor is 
also exponentially modeled except for within clouds where 
the density is set to correspond to 100% relative humidity. 
Liquid water density in clouds is taken to be triangular in 
shape similar to clouds modeled by Paris (17) . Cloud base 
height may be specified or optionally may be automatically 
determined by the point at which a parcel of air lifted 
adiabatically , with the potential temperature and mixing 
ratio constant, reaches its dew point. Cloud width and 
maximum density may also be specified manually or auto- 
matically by a quadratic relationship between cloud width 


and maximum density which was fitted to typical cumuliform 
cloud data so that specifying the columnar liquid water 
density determines the cloud size and shape. 

The atmospheric radiation, Ig^y# is determined from 
the radiative transfer equation which for a non- scattering, 
horizontally stratified atmosphere in local thermodynamic 
equilibrium is given by (2) , , - 
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(33) 


~dZ + a ^SKY = a I ATM ' 

where a is the total absorption coefficient of the atmosphere 
obtained by summing the absorption coefficients due to 

• i 

• •» 

oxygen, water vapor, and liquid water as given by equations 

| ’ 

(30), (31), and (32). I ATM i s the Plank blackbody function 
for local thermodynamic temperature T ATM , 

X ATM = (2hv 3 /c 2 )/(exp[hv/kT ATM ] - 1). (34) 

The atmospheric radiation may be specified by a brightness 

' • • : j . • • •' 

temperature T D defined by 

O 

I SKY = (2hv 3 /c 2 ) / (expIhv/kTg] - 1). (35) 

The solution of equation (33) is greatly simplified 
when the Rayleigh- Jeans approximation to the Plank blackbody 
function is used since then the intensity and brightness 
temperature of the radiation are linearly related. Both 
equations (34) and (35) are then of the form 

I = 2Kv 2 T/c 2 . (36) 

The solution of equation (33) integrated along any direction 
at an angle 0 to the vertical and for any height h in the 
atmosphere is then given by 



( 37 ) 


t b U p = I a( z ) t atm ( z) L ( h ' z ' 0 ) dZ(0) + L(h,<*>,0)T, 
* h 


and 


h 

r 


B down 


a(2)T ATM (z> L < z ' h ' 0 ) dz ( 0 > + L(O,h,0)T (38) 


Here L(a,b,0) is the transmission factor 


L(a,b,0) = exp f- I a(z)dZ(0) 


(39) 


and dz (0 ), which takes into account the earth's curvature, is 
given by 


dz (0) = 


= ! 

i - 

a_ 

sin 0 

2 



a 

+ h 

m 



- 1/2 


dz , 


(40) 


where a is the earth's radius. T c is the cosmic background 
temperature of 2.8°K (18) and T is the brightness temperature 

. . r 9 1 ' ' 


46 



of the total radiation from the ground surface. Equation 
(37) gives the radiometric brightness temperature of the 

downcoming radiation as seen looking up at an angle 0 from a 
height h in the atmosphere. Similarly equation (38) is the 
upcoming radiation .as seen looking down from height h. 

Since represents the total ground Radiation it is composed 
of both the radiation emitted by the ground and the down- 
welling sky radiation reflected from the ground. The 
downwelling sky radiation at the surface is given by equation 
(37) with h = 0. 

The Rayleigh- Jeans approximation to Plank's law is 
better than one percent for frequencies less than 100 GHz 

and temperatures greater than 250°K. However Stogryn (19) 

I ■ ■ 

has recently shown that the error introduced by using the 

j ' ' 

Rayleigh- Jeans approximation may be reduced by at least two 

orders of magnitude by changing the boundary condition 

i 

describing the downard flowing flux at the top of the earth's 
-atmosphere. This improvement in accuracy is effected by 
replacing T in equation (37) by the ficticious temperature 
T - where. 



c 2 T 

2Kv 


2 I c 


hv 
2K * 


(41) 


Here I c is the intensity of the cosmic background radiation 
and is obtained from Plank * s blackbody law (i. e . equation 



(35) with T = T = 2.8°K. This change in boundary condition 

D U 

is also reflected in T g ^ OVJn as expressed by equation (38) 

since T^ contains the reflected downwelling sky radiation at 

the surface which is given by equation (37) with h = 0. 

I SKY '*' S 9-*- ven by (35) using Tg from (37) or (38) as appropriate 

ATMOS does a Simpson's rule integration of the above 

equations to obtain either T ■ or T, at angle 9 from 
^ up down 3 

height h. The number of divisions and upper limit of the 
contributing atmosphere are set manually. 

Figure 7 is an example of the upward looking loss due 
to the separate contributing atmospheric parameters, molecular 
oxygen, water vapor and liquid water, as generated by ATMOS . 

The relative effect of the three constituents can be seen 
from 0 to 50 GHz. The calculations are for a standard 
atmosphere. The integrated columnar water vapor densities 

2 

(also called precipitable water) of 1.0, 3.5, and 6.0 gm/cm are 
typical values for a dry winter day, an average day and a 
very humid summer day, respectively, at Washington, D. C. 

2 

The liquid water columnar densities of 0.04 and 0.2 gm/cm 
represent typical stratocumulus, the majority of clouds, and 

moist cumulus clouds, respectively. As can be seen atmospheric 

• 1 • . ' 

. -I • • . ' • ... • ' / 

effects are most pronounced at frequencies above about 10 GHz. 

.. | J 

Therefore they are very important for the interpretation of 
Sr-193 data but relatively unimportant for the analysis of 
S-194 measurements. Not only are atmospheric effects small 
at the 1.4 GHz frequency of S-194 but there is very little 
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ATMOS PHRIC LOSS (PERCENT) 



FREQUENCY (GHZ) 


Figure 7 


variations from clear, dry conditions to overcast with rain. 
Only the most severe downpour results in any significant 
change at 1.4 GHz. The effective atmospheric correction for 
the S-194 measurements amounted to an average increase of 
about 5.2°K and varied only a few tenths of a degree over the 
entire range of the measurements . 

It is necessary to translate the complete gain and 
polarization properties of a satellite-borne antenna from a 
coordinate system referenced to the antenna, with respect to 
which the antenna properties are most naturally and normally 
measured, to a coordinate system oriented with respect to 
the local vertical and compass points of some observed 
location on the earth's surface, with respect to which 
radiation reflected and emitted from the earth is most 
naturally expressed. This translation will allow the power 
received by the antenna from any source of specified extent, 
intensity and polarization to be readily calculated. 

The earth-satellite geometry is shown in Figure 8. The 

A /N A 

unit vectors (i, j, k) form an orthogonal coordinate system 
with the origin at the center of the earth. The earth's 

A A 

axis lies in the (j, K) plane with north in the direction of 

1 A A 

the unit vector N. The satellite lies on the k axis at a 
height h above the earth's surface. It is moving in the 

A A 

direction of the unit vector s, which is orthogonal to k, 
at a heading H measured from north through east in the plane 
tangent to the earth at the sub-satellite point. The antenna 
is pointed in the direction of the unit vector p and views 

50 
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a point on the earth's surface located by the vector R q 
from the earth's center. The angle of incidence on the 
earth's surface is given by I Q . The variable point for 
integration over the antenna beam is located at R in the 

r 

A 

direction u with an incidence angle on the surface of I . 

P P 

The latitude and longitude of the sub-satellite point, the 

antenna look point, and the integration point are given by 

(L , l ) , (L , l ) and (L , 1 ), respectively, 
s s o o P P 

The satellite is located at R g and letting the radius 
of the earth be given by a, we have, 



A 


S 

A 

N 


A 


(a + h) k, 

A A 

sin Hit cos H j 

A 

cos L j + sin L 
s J s 



(42) 

(43) 

(44) 


R q = R s - Ua+h) y Q -k + ^a 2 - (a+h) 2 [1- (y Q . k) 2 ] }y Q , (45) 
cos I Q = - y Q • R Q /a, (46) 

A 

sin L q = N • R Q /a , (47) 

A A 

cos (.5, -£. ) = i • (N x R o /a)/cos L q . (48) 


-i ... ^ 

Expressions for R^, 1^, and 2-^ are given by changing the 
subscripts o to p in equations (45), (46), (47) and (48). 

The antenna properties may be conveniently specified in 
a spherical coordinate system referenced to the antenna and 
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expressed in terms of Stokes parameters as shown in Figure 9. 

A /N A 

The unit vectors (y o# y^, y v ) form an orthogonal coordinate 

A A /N 

system as do the unit vectors (y , y 0 , y^)„ The direction 
of maximum response of the antenna would generally be in the 

A 

direction of y with the aperture plane coinciding with the 

/V A 

y y ) plane. The Stokes parameters are given in normalized 
form as 

1^(6, <j>) = f(6,<M 

a 

Q_ (6,4)) « f ( 6 , $ ) cos (23) cos {2yJ (49) 

a 

U (0,<}>) = f ( e , 4? ) cos (23) sin (2x) 

a ■ 

V ( 6 , <}>.) = f (0,$) sin (23) , 

* ® 

where f(o,o) = 1, x specifies the orientation of the major 
axis of the polarization ellipse and tan 3 is equal to the 
ratio of the minor axis to the major axis. A "perfect" 
linearly polarized antenna would have 3=0. It might be 
described to be "vertically" polarized when the electric 
vector of its radiation pattern coincides with the meridian 
planes of a sphere with a polar axis in the direction y . 

In this case cot x = cos 9 tan 4>. Naturally I (e,4>) 
must be zero in the direction of the poles since the meridians 
coverge at the poles and x v is not uniquely given there. A 
"horizontally" polarized antenna would be specified by 
having the electric vector of its radiation pattern ortho- 
gonal to the meridian planes; i.e. = X y + tt/2. 
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The Stokes parameters referenced to a second coordinate 

A A A /\ 

system, indicated by primes, where p' = p but pA and p! 

Jr u t 

1 A 

aire rotated, in the positive sense, about p by an angle Ax 
! P 

A A 

from p n and p, , are related to those of the unprimed system 
u 9 

by 


I’ = I 

Qj ' = Q cos (2Ax) + U sin (2Ax) 


(50) 


U 1 ’ = U cos (2Ax) ~ Q sin (2Ax) 


V' = V. 


The rotation Ax necessary to reference the antenna Stokes 

i a 

parameters to a coordinate system where p^ is perpendicular 

A 

to, and p^ lies in, the plane of incidence at the integration 
point on the earth's surface is given by 


cos Ax = P fi * [R x p ]/Ia sin I ] 


(51) 


AAA 


The unit vectors (p , p„, p.) are given in terms of 

p o <p 


A A A 


<>V V M v> by ' 


y = sin 6 cos <f> p h + sin 9 sin 4> p v + cos 0 p ( 


y_ = cos 9 cos d p, + cos 9 sin <p p„ - sin 0 p^ (52) 
^6 n v o 


y^ = - sin <f> y h + cos 4> y y . 
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It is now necessary to relate the antenna coordinate 
system to the coordinate system fixed in the satellite. 
This may be done by the use of modified Euler’s angles (4^, 
\l>2' ^3) as shown in Figure 10 . Here the line of Nodes is 

A A 

given by k x and 


A 

[ sin i>2 sin 1 p^ - cos 4 ']_ cos 4^ cos ] s 

A A 

+ [cos ^2 sin 4*3 + cos sin ip 2 cos ^3] k x s 

A 

+ [- sin ip^ cos 4^3 k 


+ [cos 4 j 2 cos ip ^ ~ cos 1 p^ 

A 

+ sin 4^ sin 4^] k 


P Q = [sin 4>2_ cos 4^3 s 

A /N 

+ [- sin ipj sin 4 i 2 3 k x s 

A 

+ [-cos 4 ^] k. ( 55 ) 

’ . ' . ‘ ••• • - ..... i-. . : • 

Note that if the antenna is scannned along the in-track 

A A A 

direction, that is y Q is scanned in the (s , k) plane, 

4 ; x = INTRACK SCAN ANGLE 

ir/2. 


^2 ” 
% = 
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Figure 10 
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For cross-track scanning, where y Q is scanned in the 

A A A 

(k, k x s) plane, 

^1 = ^c = CR0SSTRACK SCAN angle 

= " tt/2 (57) 



For scanning along the surface of a cone with its axis in the 

A 

-k direction, 

= APEX ANGLE OF CONE 

\p 2 ’= CONE SCAN ANGLE (58) 

^2 = Tr /2 • 

W • 

Substituting (43), (53), (54), and (55) into (52) gives y in 

. ! i !" ‘Jr.. 

the earth centered coordinate system and expressed in terms 
of Euler's angles and the antenna integration angles, 


= [cosSsinip^sin (Ht^) - sinBcos (H+i/^) sin (<f> + ^ y) 

. ■ ;• i : ■ ’ . ... : - a 

- sin0costj^sin (H+i^) cos ( ) 3 i 

+ [cos0sin^cos (H+ij^) + sin0sin (H+iJ;^) sin (59) 

A 

- sin0cost^cos (H+iJi^) cos ((fi+ij;^) ] j 

A ' .. 

+ [- cosQcos^ - sin0sin^^cos ($+^) j k. 
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Hg may be obtained by replacing 0 with (0+tt/2) in (59) 

A /\ 

and u becomes u when 0=0. 
p o 

Given the orientation of the antenna with respect to 

the satellite \p 2 , and ’j> 3 , the satellite coordinates, 

altitude and heading, L , Z , h and H and the running variables 

of integration 6 and <f> , the angle of incidence and coordinates 

at the look point I , L , and l and the angle of incidence 

o o o 

and rotation angle for the Stokes parameters at the integra- 

" . j ■ • • ' ■; . • • - ... . ' 

tion point I and Ax may be found from equations (42) , (44) , 

r 

(45), (46), (47), (48), (51) and (59) to be: 


sin I 

o 

sin L 

o 


(1 + h/a) sin (60) 

sinL cos (I -h. ) + cosL sin (I -ip-,) cos(H+bo) (61) 
s o • 1 s o 1 Y 2 


cos (1 -l ) = [cosL 

O U b 


sin Ip = (1 + h/a) 


COS (Ig,"^) 

- sinL g sin(I Q -i|;^) 


- [cos0cosip^ + 


(62) 

cos (H+iJ^ ) ] /cosL q , 


2 

sinOsimjj^cos (cj)+^ 3 ) ] (63) 


cos Ax= (1 + h/a) simp, sin((p+b 3 )/sin I . (64) 

The Stokes parameters of the radiation incident on the 
antenna are found by summing the Stokes parameters of the 
radiation emitted and reflected from the sea and emitted by 
the atmosphere and applying the appropriate atmospheric 
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attenuation. Let the Stokes parameters of the total radiation 
be given by [I t , Q t r U t , V^] and the Stokes parameters of 

i 

the antenna, properly translated, be represented by [I , Q , 

ci cl 

U . VJ, then the power received by the antenna per unit 

Cl ci 

frequency interval is, e.g. (20) , 

P - <V 2)/ ‘Va + Q t Q a + U t U a + V t V a f d!2 ' (65) 

where A is the effective area of the antenna. The effective 
e 

area may be found from the normalized antenna response 

pattern f(0,4>) and a measurement of the radiation efficiency 

n of the antenna or a measurement of the maximum gain of the 

antenna G where 
o 

«r. i I 

A e - nl 2 //f = A 2 G 0 /4 tt . (66) 

Measurement of n is generally very difficult in practice 
and, since it is usually very close to unity, it is often 
taken to be one. 

The power received by the antenna is often specified in 
terms of the so called antenna temperature T^. This is the 
temperature to which an impedance, equal to the impedance of 
the antenna, would have to be raised in order that when it 
is substituted for the antenna it deliver the same power to 
the input of the receiver, as does the antenna. The antenna 
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temperature is related to the power received per unit 
frequency interval by 


T a = P/k 


( 67 ) 


Equation (65) is then the desired integration of the 
total radiation incident on the antenna over the antenna 
reception pattern. It represents the principal output of 
ANTTEMP and is quite general. 

Solar radiation reflected from the sea surface can be 
appreciable; especially at 1.4 GHz. It is strongly dependent 
on the ocean surface roughness. In order to calculate the 
."sun glint" scattered into the antenna from the ocean surface, 
the rough surface scattering or geometries optics model as 
developed by Stogryn (10) was used. Following Stogryn, and 
assuming the radiation from the sur. is unpolarized with a 
brightness temperature the brightness temperature of 

the horizontally polarized radiation scattered in the direction 
(0,4>) is given by 


T h 4 if ^ T sun v, hh ’ 'hv J 


(Yv,v> + YwJ dfl • 


( 68 ) 


Since (Y hh + Y hv ) changes very little over the relatively very 
small solid angle subtended by the sun 


, sun fisun . . . 

fc h 4? lY hh + T hv> ' 


( 69 ) 
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where ft is the solid angle of the sun. The scattering 
sun 


coefficient is given by 


2 2 
(a +r> 


hh 


hv 


<r 2 cos 2 A + r 2 s in 2 A) e 2 9 V 
2g B cos 0 n 


( 70 ) 


2 2 

where R. and R are the Fresnel reflection coefficients 
h v 

evaluated at an angle i and 


cos i = j (ct 2 + 3 2 + B 2 ) 1 ^ 2 , 


a = sin 0 cos <J> + sin 0„ . , 

sun 


3 = sin 6 sin <J> , 


B = cos 0 + cos 0 , 

sun ' 


. 2 . 2 nr 

« sm 6 sm 0^ 

_ . 2, sun 

Sm A = - 


sin 2 (2i) 


g = 0.0015 + 0.00131 W. 


(71) 


Here 9 sun is the incidence angle of the sun on the mean surface 
and W is the marine wind speed in knots. The mean-square surface 
slope, g , was obtained from Cox and Munk (21), (22) by averaging 

the mean-square slopes in the upwind-downwind and crosswind 
directions that they derived from measurements of the sun's 
glitter pattern on the ocean. The brightness temperature of 
the vertically polarized radiation, T , is obtained from 
equations (69) and (70) by interchanging the subscripts h and v 
throughout. The Stokes parameters for the scattered solar 
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radiation are given by, 


(72) 

Adopting one-half degree for the diameter of the solar disc and 

'5 ' i 

10 °K for the solar brightness temperature at 1.4 GHz (20), the 

antenna temperature increase due to sun glint was calculated 
using ANTTEMP and is given in Figure 11 as a function of the 

t ' ’ - ' 

solar incidence angle on the surface for several surface wind 

speeds. As can be seen, the effects of sun glint decrease 

rapidly with solar incidence angle and wind speed but can be 

significant. These calculations, though serving to show the 

j ' .... ' 

qualitative dependence on solar incidence angle and wind speed, 

can only be regarded as approximate because of the uncertainty 

in the scattering coefficient. There is some evidence (8) 

that the effective mean- square surface slope at 1.4 GHz is only 

one- third of that given by equation (71) . This results from 

surface structure smaller than about one observational wavelength, 

22 cm at 1.4 GHz, being ineffective in scattering 1.4 GHz 

radiation. Since the short surface wavelengths contribute 


= 


K 


2 (T h + V 


Q s - fr (T h - V 


u = 0 , 


V_ = 0 . 
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disportionately to the mean-square surface slope, its effective 
v^lue decreases with decreasing observational frequency. The 
reduction at 1.4 GHz to one-third of the "optical" value is 
similar to the reduction measured by Cox and Munk for an oil 
slick covered surface compared to a clean surface. Since they 
calculated that an oil slick damped surface waves of 30 cm 
wavelength by a factor of 10 and that waves shorter than this 
were essentially eliminated, the sea surface at 1.4 GHz may 
very well be represented by the mean-square slope of a slick 
covered surface. Using a mean-square slope one-third of that 
given in equation (71) results in antenna temperature increases 
due to sun glint two to three times larger than the values given 
in Figure 11. 

ANALYSIS OF S-194 DATA: 

A compendium of all the S-194 measurements analyzed is 
presented in Table 2. This consists of 12 0 ind'ependent 
observations over 27 EREP passes and includes data over the 
Equatorial Pacific, North Pacific, North Atlantic, Equatorial 
Atlantic, and South Atlantic Oceans, the Gulf of Mexico, 
Caribbean Sea, the Great Lakes Region, the Southwest American 
Desert, and the Amazon River Panalopy. In general, over the 
open ocean regions a measured antenna temperature consists 
of an average over a statistically stable interval lasting 
for a period of several minutes. When the antenna response 
is varying rapidly as in Pass 14 over the Great Lakes and in 
regions of land-water boundaries, precise positions and 
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DATE 

SKYLAB 
PASS # 

GMT . .. 

LAT. . 
°N 

LONG. 

°W 

SUN 

ELEV 

ANG 

(°) 

6/5/73 

5 

18:04: 00 

24.5 

92.0 

88 

6/10/73 

7 

14:29:00 

32.0 

76.6 

54 



14:35:00 

15.4 

60.0 

69 



14:40:00 

0.9 

48.8 

87 

6/11/73 

8 

15:19:00 

33.0 

96.8 

48 



15:22:00 

25.9 

88.. 5 

55 



15:24:00 

10.9 

82.0 

61 



15:31:00 

-1.3 

66.0 

62 

6/12/73 

9 

13:05:00 

29.2 

63.8 

47 

• 


13:09:00 

15.6 

50.8 

58 



13:10: 00 

13.9 

49.3 

59 



13:13:00 

6.4 

43.5 

61 



13:15:00 

-0.2 

38.5 

60 



13:16:00 

-2.4 

36.8 

59 



13:18:00 

-9.8 

31.2 

55 



13:19:00 

-11.8 

29.6 

54 

6/13/73 

10 

14:01:00 

12.6 

67.3 

53 

8/4/73 

13 

17:09:00 

49.7 

160.6 

16 


TABLE 2 



PPT 

°C 

KTS 

«°K) 

(°K) 

( °K) 


36 

27 

12 

114.2 

103.7 

-.- 

Sun Glint 

36 

26 

5 

94.7 

95.0 

-0.3 

Atlantic 

36 

27 

12 

101.1 

100.0 

— 

Sun Glint 

21 

29 

4 

118.8 

118.6 

— 

Sun Glint 

— 

— 

— 

251.9 

— 

— 

South West Desert 

36 

28 

3 

93.1 

94.0 

-0.9 

Gulf of Mexico 

36 

29 

8 

94.7 

95.0 

-0.3 

Caribbean 

— 

— 

— 

271.6 

— 

— 

Amazon Jungle 

36 

27 

5* 

93.7 

94.8 

-1.1 

North Atlantic 

36 

27 

10* 

95.3 

95.9 

-0.6 

North Atlantic 

36 

27 

10* 

96.2 

95.9 

.3 

North Atlantic 

35.5 

• 27 

10* 

96.9 

96.2 

.7 

North Atlantic 

36 

28 

5* 

95.9 

94.6 

1.3 

Equatorial Atlantic 

36 

28 

5* 

96.5 

94.6 

2.3 

South Atlantic 

36.5 

27 

10* 

95.6 

95.3 

.3 

South Atlantic 

36.5 

27 

5 

94.7 

94.5 

.2 

South Atlantic 


35.5 27 12 94.6 96.5 -1.9 


33 


12 


5* 


95.9 97.6 -1.7 


Caribbean 
North Pacific 




8/4/73 

13 

17:09:00 

50.2 

149.2 

24 


■ ; : 1 T ; 

17:11:00 

49.4 

137.8 

32 



17:13:00 

33. 0 

96.0 

68 



17:22:00 

25.4 

87.3 

79 



17:25:00 

-1.3 

65.2 

65 

8/5/73 

14 

15:01:00.45 

44.891 

89.470 

43 



15:01:15.46 

44.468 

88.393 

44 



15:01:24.09 

44.200 

87.662 

44 



15:01:32.35 

43.949 

87.036 

45 



15:01:40.96 

43.694 

86.389 

45 



15:01:49.56 

43.486 

85.744 

46 



15:02:05.1 

42.955 

84.614 

47 



15:02:24.95 

42.334 

83.186 

48 



15:02:33.22 

42.068 

82.602 

49 



15:02:41.16 

41.812 

82.059 

49 



15:02:50.09 

41.512 

81.426 

50 



15:02:58.69 

41.224 

80.835 

50 



15:03:60.63 

40.955 

80.291 

51 



15:03:37.39 

39.887 

78.237 

53 



15:04:08.15 

38.774 

76.248 

55 



15:04:16.43 

38.471 

75.702 

55 



15:04:25.02 

38.151 

75.183 

56 



15:04:33.96 

37.816 

74.638 

56 


2 (continued) 


MEAS. CALC. AT 


L.INITY SST 
?PT °C 


SPEED TEMP 
KTS ( °K) 


ANT. 

TEMP 

(°K) 


MEAS.' 

CALC. 

(°K) 


COMMENTS 


33 13 


97.5 -1.9 


North Pacific 


32 

] 5 

5* 

96.7 

98.7 

-1.5 

North Pacific 

— 

— 

— 

270.2 

— 

— 

South West Desert 

36 

28 

9 

109.2 

102.9 

— 

Sun Glint 

— 

— 

— 

274.5 


— 

Amazon Jungle 

— . 

— 

— 

265.6 

264.8 

0.8 

• Wisconsin 

— 


— 

246.6 

243.8 

2.8 

Wisconsin 

— 

— 

— 

195.0 

192.5 

2.5 

Wisconsin 

0 

12 

10 

153.0 

155.3 

-2.3 

Lake Michigan 


— : 

— 

189.2 

191.5 

-2.3 

Michigan 

— 

— 

— 

242.7 

244.0 

-1.3 

Michigan 

— 

— 

— 

259.9 

262.1 

-2,2 

Michigan 

— 

— 

— 

229.4 

232.0 

-2.6 

Michigan 

— 

— 

— 

190.2 

191.2 

-1.0 

Michigan 

0 

23 

10 

173.0 

172,0 

1.0 

Lake Erie 


— 


208.4 

207.5 

0.9 

Ohio 

— 

— 

— 

252.2 

250.9 

1.3 

Ohio 

— 

— 

— 

262 . 6 

263.6 

-0.9 

Ohio 

— 

— 

— 

266.5 

266.1 

0.4 

Pennsylvania 

— 

— 

— 

223.6 

221.1 

2.5 

Maryland 

— 

— 

■ — 

203.5 

200.8 

2.7 

Maryland 

“ 

— 

— 

153.9 

153.9 

0 

Maryland 

— 

— 

— 

110.0 

110.8 

-0.8 

Coast of Maryland 



TABLE 2 (continued) 








SUN 




MEAS . 

CALC. 

AT 



SKYLAB 





ELEV. 




ANT . 

ANT. 

MEAS.- 


DATE 

PASS // 

GMT 


LAT. 

LONG. 

ANG. 

SALINITY 

SST 

SPEED 

TEMP 

TEMP 

CALC. 

COMMENTS 





°N 

°W 

(°) 

PPT 

°C 

KTS 

(°K) 

(°K) 

(°K) 

• ; ; 

8/05/73 

14 

15:04:42. 

84 

37.478 

74.078 

57 

35.5 

26 

4 

97.6 

96.3 

1.3 

North Atlantic 



15:05:01. 

73 

36.253 

72.172 

59 

36.0 

27 

5 

94.4 

94.5 

-0.1 

North Atlantic 

8/05/73 

15 

16:41:00 


28.9 

86.2 

70 

35.0 

29 

5* 

104 . 2 

98.8 

— 

Suti Glint 



16:46:00 


15.0 

73.0 

87 

36.0 

28 

5 

104.4 

107.8 

— 

Sun Glint 

8/08/73 

16 

15:55:00 


59.5 

143.4 

15 

32.0 

14 

5* 

97.4 

98.1 

-0.7 

North Pacific 



15:57:00 


47.9 

133.3 

22 

32.0 

16 

5* 

98.3 

98. 2 

0.1 

North Pacific 



16:05:00 


32.8 

100.3 

50 

— 



283.3 

— 

. 

Southwest Desert 



16:08:00 


24.2 

91.5 

59 

36.0 

29 

5 

94.3 

94.2 

0.1 

Gulf of Mexico 

8/09/73 

17 

13:51:00 


26.8 

64.8 

53 

36.5 

28 

10* 

92.4 

93.7 

-1.3 

North Atlantic 



13:59:00 


4.7 

45.8 

70 

36.0 

27 

5* 

93.9 

94.8 

— 

Sun Glint 

9/02/73 

23 

17:55:30 


10.7 

58.6 

60 

32.0 

28 

5* 

100.0 

99.8 

0.2 

Off Northeast Coast 
of Venezuela 



17:55:40 


11.1 

58.4 

59 

32.0 

28 

5* 

99.0 

98.0 

1.0 

North Atlantic 



17:56:20 


13.1 

56.7 

58 

35.0 

28 

20 

98.1 

97.9 

0.2 

North Atlantic 

rn 

is! s 


17:57:10 


15.7 

54.6 

55 

36.0 

28 

30 

99.0 

98.8 

0.2 

Closest approach to 
tropical storm 
Christine 

t-p* W. 

Ps 


17:57:20 


16.3 

54.1 

54 

36.0 

28 

25 

98.1 

98.0 

0.1 

North Atlantic 


17:58:00 


18.0 

52.6 

53 

36.0 

28 

10 

95.9 

95.6 

0.3 

North Atlantic 

; • 

1*^*4 

: Cv *~g 


17:59:00 


21.1 

49.9 

49 

36.0 

28 

5 

94.9 

94.8 

0.1 

North Atlantic 


18:00:00 


23.4 

47.6 

47 

37.0 

28 

5 

94.0 

94.1 

-0.1 

North Atlantic 















o , 














§5J 

SI 















TABLE 2 (continued) 


DATE 

SKYLAB 
PASS // 

GMT 

LAT. 

*N 

LONG. 

°W 

SUN 

. ELEV. 
ANG. 
(°) 

SALINITY 

PPT 

SST 

°C 

SPEED 

KTS 

MEAS. 
ANT. 
TEMP 
( ° K) 

CALC. 

ANT. 

TEMP 

(°K) 

AT 

MEAS.- 

CALC. 

(°K) 

COMMENTS 

9/09/73 

29 

19:09:00 

12.8 

118.0 

79 

33.0 

25 

0 

115.4 

111.2 

— 

Sun Glint 



19:27:00 

50.1 

48.0 

20 

32.0 

11 

14 

98.4 

99.6 

-1.2 

North Atlantic 



19:30:00 

49.4 

30.7 

9 

35.5 

16 

11 

98.1 

97.8 

0.3 

North Atlantic 

9/10/73 

31 

20:05:00 

29.4 

121.2 

64 

34.0 

20 

5* 

98.0 

97.5 

0.5 

North Pactific 

9/13/73 

38 ' 

18:09:00 

50.2 

48.8 

30 

33.0 

10 

10* 

99.9 

98.6 

1.3 

North Atlantic 



18:10:00 

50.1 

44.0 

28 

34.5 

13 

10* 

100.8 

98.1 

2.7 

North Atlantic 



18:12:00 

49.1 

32.8 

21 

35.0 

17 

10* 

100.1 

97.9 

2.2 

North Atlantic 

9/14/73 

40 

17:05:00 

5.4 

123.3 

44 

35.0 

27 

5* 

99.5 

95.4. 

-0.9 

Equatorial Pacific 



17:08:00 

14.2 

116.3 

51 

34.0 

28 

5* 

94.7 

95.9 

-1.2 

North Pacific 



17:09:00 

18.9 

112.4 

53 

34.0 

28 

5* 

93.8 

95.9 

-2.1 

North Pacific 



17:30:00 

48.0 

21.4 

21 

35.5 

16 

10* 

94.3 

97.7 

-3.4 

North Atlantic 

9/15/73 

43 

18:01:00 

22.2 

128.2 

50 

34.0 

23 

6 

98.3 

96.8 

1.5 

North Pacific 



18:02:00 

25.5 

125.0 

51 

34.0 

22 

6 

98.3 

96.9 

1.4 

North Pacific 


.L 

18:04:00 

30.4 

119.7 

52 

• 

34.0 

19 

8 

100.2 

97.2 

3.0 

North pacific 

12/04/73 

60 

16:45:00 

49.3 

150.3 

-9 

33.0 

6 

5* 

95.3 

97.4 

-1.® 

North Pacific 

1/04/74 

73 

19:40:00 

47.7 

13.7 

23 

35.5 

12 

10* 

96.2 

96.6 

-0.4 

North Atlantic 

1/06/74 

74 

17:57:00 

27.2 

88.1 

40 

36.0 

23 

5* 

94.5 

95.6 

-1.1 

Gulf of Mexico 



17:58:00 

28.5 

86.7 

39 

36.0 

21 

5* 

95.5 

95.9 

-0.4 

Gulf of Mexico 



18:07:00 

47.9 

48.6 

10 

34.0 

2 

. 10* 

96.4 

97.6 

-1.2 

North Atlantic 



18:11:00 

50.2 

26.2 

-3 

34.0 

2 

10* 

97.0 

97.6 

-0.6 

North Atlantic 


TABLE 2, (continued) 


DATE 

SKYLAB 
PASS f 

GMT 

LAT. 

°N 

LONG. 

°W 

SUN 

ELEV. 

ANG. 

O 

SALINITY 

PPT 

SST 

°C 

SPEED 

RTS 

MEAS. 

ant. 

TEMP 

(°K) 

CALC. 

ANT. 

TEMP 

(°R) 

AT 

MEAS.- 

CALC. 

(°R) 

COMMENTS 

1/06/74 

74 

18:13:00 

49.5 

14.8 

-9 

35.5 

11 

10* 

96.1 

97.7 

-1.6 

North Atlantic 

1/07/74 

76 

17:14:00 

25.4 

85.3 

42 

36.0 

23 

10* 

96.9 

96.7 

0.2 

Gulf of Mexico 



17:24:00 

46.7 

48.7 

16 

33.0 

1 

10* 

97.6 

97.8 

-0.5 

North Atlantic 

1/08/74 

78 

16:29:00 

19.8 

85.7 

43 

36.0 

23 

5* 

95.5 

95.6 

-0.1 

Gulf of Mexico 



16:32:00 

27.8 

77.9 

38 

36.5 

24 

4 

94.2 

95.3 

-0.9 

North Atlantic 



16:34:00 

31.7 

73.4 

36 

36.5 

24 

4 

95.8 

95.3 

0.5 

North Atlantic 



16:35:00 

34.9 

69.1 

33 

36.5 

20 

4 

95.8 

95.9 

-0.1 

North Atlantic 



16:36:00 

37.3 

65.6 

30 

36.5 

20 

10 

96.5 

96.9 

-0.4 

North Atlantic 



16:37:00 

39.2 

62.5 

28 

36.0 

17 

10* 

98.4 

97.4 

1.0 

North Atlantic 



16:38:00 

40.6 

59.8 

27 

36.0 

13 

10* 

97.1 

97.5 

-0.4 

North Atlantic 



16:39:00 

43.2 

54. 2 

23 

34.0 

3 

10* 

97.0 

97.7 

-0.7 

North Atlantic 



16:40:00 

44.4 

51.3 

22 

33.0 

3 

7 

96.1 

97.4 

-1.3 

North Atlantic 



16:41:00 

47.6 

40.4 

18 

35.5 

12 

10* 

94.8 

97.7 

-2.9 

North Atlantic 

1/09/74 

7? 

15:51:00 

31.1 

69.4 

35 

* 36.5 

23 

8 

96.4 

96.1 

0.3 

North Atlantic 



15:52:00 

33.8 

62.1 

33 

36.5 

20 

6 

97.3 

96.6 

0.7 

North Atlantic 


/ 

15:55:00 

29.6 

56.9 

28 

36.0 

14 

30 

101.5 

100.7 

0.8 

North Atlantic 



15:57:00 

44.4 

46.4 

23 

36.0 

6 

48 

103.0 

103.3 

-0.3 

North Atlantic 



15:58:00 

46.3 

40.9 

20 

36.0 

12 

48 

103.9 

102.5 

1.4 

North Atlantic 



15:59:00 

47.2 

37.5 

18 

36.0 

12 

48 

103.6 

102.5 

1.1 

North Atlantic 



16:00:00 

48.7 

30.4 

15 

35.5 

10 

20 

100.5 

99.3 

1.2 

North Atlantic 


REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 


TABLE 2 (continued) 


DATE 

SKYLAB 
PASS if 

GMT 

LAT. 

°N 

LONG. 

°W 

SUN 

ELEV. 

ANG. 

(°) 

SALINITY 

PPT 

SST 

°C 

SPEED 

KTS 

MEAS. 

ANT. 

TEMP 

(°K) 

CALC. 
ANT. ' 
TEMP 
(°K) 

AT 

MEAS.- 

CALC. 

(°K) 

COMMENTS — 

1/09/74 

79 

16:02:00 

50.1 

18.4 

9 

35.5 

10 

10* 

99.2 

97.7 

1.5 

North Atlantic 



16:04:00 

50.2 

11.0 

6 

35.5 

11 

13 

97.9 

97.8 

0.1 

North Atlantic 

1/11/74 

81 

17:49:00 

46.8 

24.2 

6 

36.0 

12 

10* 

97.7 

97.5 

0.2 

North Atlantic 

1/27/74 

92 

12:18:00 

13.3 

143.4 

53 

34.0 

26 

14 

99.3 

98.7 

0.6 

North Pacific 



12:20:00 

19.5 

138.1 

49 

34.0 

22 

20 

100.3 

100.4 

-0.1 

North Pacific 



12:22:00 

23.3 

134.6 

46 

34.0 

21 

23 

99.3 

100.1 

-0.8 

North Pacific 



12:23:00 

27.4 

130.4 

42 

34.0 

18 

17 

99.0 

99.4 

-0.4 

North Pacific 



12:43:00 

44.7 

36.6 

22 

36.0 

14 

45 

104.3 

103.0 

1.3 

North Atlantic 



12:46:00 

38.5 

23.4 

31 

36.0 

17 

31 

100.5 

100.6 

-0.1 

North Atlantic 



12:47:00 

35.8 

19.1 

35 

36.0 

17 

16 

99.2 

98.3 

0.9 

North Atlantic 

1/29/74 

95 

17:39:00 

24.0 

90.4 

47 

36.0 

25 

11 

97.9 

96.5 

1»4 

North Atlantic 



17:49:00 

-3.3 

68.3 

69 

— 

— 

— 

270.3 

— 

— 

Amazon Jungle 

1/30/74 

96 

16:44:00 

49.6 

139.0 

-2 

33.0 

6 

17 

98.3 

99.3 

-1.0 

North Pacific 



16:45 : 00 

48. 7 

132.4 

3 

33.0 

7 

11 

99.1 

98.4 

0.7 

North Pacific 


w. ' ..... . 

16:54:00 

33.1 

96.0 

34 

— 

— 

— 

198.3 

- — 


Southwest Desert 



16:57:00 

24.7 

86.3 

45 

36.0 

26 

5 

98.6 

95.2 

3.5 

North Atlantic 



16:59:00 

16.8 

79.1 

55 

36.0 

26 

18 

99.5 

97.4 

2.1 

North Atlantic 



17:06:00 

-1.3 

65.0 

72 



— 

275.1 

■- — 

— 

Amazon Jungle 


TABLE 2 (continued) 


SUN “ MEAS. CALC. AT 

SKYLAB ELEV. ANT. ANT. MEAS.- 


DATE 

PASS it 

GMT 

LAT. 

°N 

LONG. 
; °W 

ANG. 

~ C) 

SALINITY 

PPT 

SST 

°C 

SPEED 

KTS 

TEMP 

(°K) 

TEMP 

(°K) 

CALC. 

(°K) 

COMMENTS 

2/01/74 

98 

16:51:00 

50.1 

161.6 

-14 

33.0 

5 

5* 

98.5 

97.1 

1.4 

North Pacific 



16:53:00 

48.7 

146.7 

-4 

33.0 

7 

25 

99.2 

100.1 

-1.6 

North Pacific 



16:56:00 

45.8 

134.8 

5 

33. 0 

9 

20 

99.2 

100.1 

-0.9 

North Pacific 


* Denotes ground truth data taken more than five hours away form measurement time. 


times are recorded and averages of only a few seconds duration 
are used. 

The first column of Table 2 gives the date that the 
measurement was taken. The second column gives the Sky lab 
pass number using the continuous numbering system beginning 
with Skylab II rather than the alternative system, occasionally 
used elsewhere, which starts anew at the beginning of each 
Skylab mission. The sun elevation angle is the complement 
of the angle of incidence of the sun's rays at the point 
directly below the S-194 antenna. 

Salinity measurements in parts per thousand were obtained 
from the hundred year averages compiled by the U.S. Naval 
.Oceanographic Office (23) , and as given by Sverdrup', Johnson 

and Fleming (24). These values, while accurate on the 

| 

average, may not apply in specif ic- cases . For example, a 
recent heavy rain can result in significant temporary local 
variations in the surface salinity. It is not likely that 
this would result in an appreciable effect at 1.4 GHz however, 
except under unusually heavy rain conditions, since the skin 
depth at 25°C sea temperature increases from 0.9 cm for 35 PPT 
salinity to over 10 cm for fresh water. Therefore, S-194 
would respond to water well below the thin fresh water 
surface. Turbulance and mixing would, of course, tend to 
prevent a significant salinity gradient with depth. In 
addition, the footprint of the S-194 on the earth's surface 
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has a 60 nautical mile diameter, so that it would be necessary 
for the heavy rainfall to occur over a large area for its 
effect on the antenna temperature to be noticeable. 

Sea temperature measurements were obtained from the 

Fleet Numerical Weather Central (FNWC) , Monterey, California, 

l ......... 

except for the Great Lakes measurements which were provided 
for Pass 14 by the Detroit Office of NOAA. Wind speed in 
knots was obtained from the Daily Weather Maps produced by 
NOAA (e.g. 25) for the Great Lakes. Over the ocean wind 
speed values came from FNWC. Wind speed value labeled with 
an asterisk are estimates using data taken more than five 
hours away from the measurement time of the antenna temperature 
. The measured antenna temperature recorded in Table 2 
incorporates the corrections applied to the electronics box 
constants C 2 and C^ as discussed previously. The calculated 
antenna temperature is the output of program ANTTEMP and 
incorporates the specification of sea temperature, sea 
salinity, wind speed, and the atmospheric conditions obtained 
from all available ground truth. Although antenna temperatures 
were computed for data points with sun elevation angles in 
excess of 65° (incidence angles less than 25°), the measured 
antenna temperature minus the calculated antenna temperature 
(AT) was not calculated and the data were not used in any 
further analysis due to the extreme sensitivity of the 
antenna temperature to wind speed at these low solar incidence 
angles as shown in Figure 11. The "comment" column identifies 
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sun glint problems and gives the general location of each 
measurement. 

The first objective in the analysis of S-194 data is to 
examine those measurements taken over clear, calm ocean 
conditions to confirm the reduction and analysis techniques 
under the most ideal ocean environment before proceeding to 
more complex conditions. Three areas with calm seas, defined 
as wind speed less than or equal to five knots, minimum 
atmospheric loss dnd with very good ground truth available 
were selected. One area is in the South Atlantic and was 
measured during Pass 9 on 6/12/73. The other two areas are 
in the North Atlantic and were measured during Pass 23 on 
9/2/73. The agreement between measured and calculated 
antenna temperatures was excellent as summarized below: 


DATE 

SL PASS 

MEASURED 

CALCULATED 

. t 

i — ' 

ANT. TEMP. 

ANT. .TEMP . 

6/12/73 

9 

94.7 

94. '5 

9/2/73 

23 

94.9 

94.8 

9/2/73 

23 

94.0 

94.1 


We feel that this confirms the accuracy of the calibration 
and analysis techniques and establishes the basis for the 
investigation of more complex environments . 

The second objective in the analysis is to devise a 
method for computing the antenna temperature when the beam 
includes a land-water interface or other such areas of 
rapidly changing signal. The half-power beam-width of the 
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S-194 antenna is 15°. As can be seen in Figure 3, the main- 
beam, defined as the beam included within the first null at 
20 ° off the main-beam axis, contains 98% of the antenna 
solid angle. Further, 99% of the power received by S-194 is 
included within a 41° angle from nadir. A program MAPS was 
written to map the S-194 beam pattern out to 41° off nadir 
onto the earth's surface. MAPS approximates the beam pattern 
by a grid of 1728 points. The height above the earth's 
surface of Skylab and the latitude and longitude of the 
nadir point are inputs to MAPS. Longitude and latitude is 
computed every ten degrees in azimuthal angle and every half 
degree angle off nadir out to eighteen degrees (97.5% power 

included) then every two degrees out to forty-one degrees. 

* 

From the output of MAPS superimposed on a world map each 
point was then identified as either being on land or water 
and inputs to ANTTEMP adjusted accordingly. 

The first such model was made of an arbitrary straight 
line interface between land at a brightness temperature of 
2;76.5°K and the ocean at a brightness temperature of 96.2°K 
to examine the effects of worst case pointing error on 
antenna temperature. A study of the Skylab attitude pointing 

error which used data from the Skylab camera systems provided 

"i • ■ . . 

3a uncertainties of 2.5° for EREP Passes 1 through 20 and 
0.8° for all subsequent EREP passes. These uncertainties 
correspond to a worst case error of the center of the 
antenna footprint on the earth's surface of 10.2 nautical 
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miles for EREP Passes 1-20 and 3.2 nautical miles for subsequent 
passes. If it is assumed that the nadir look angle is 
exactly on the interface, a pointing error of 2.5° will 
result in an error in measured antenna temperature of 27.4° 

(a 14.7% error); and a pointing error of 0.8° will result in 
an error in the measured antenna temperature of 7.1°K (a 
3.8% error). It is thus clear that pointing error can be a 
significant source of error in measured antenna temperature 
if the main beam of S-194 includes such a land-water interface. 


Such an interface occurred frequently during Pass 14, 
Ground Track 61/62, August 5, 1973. This Pass started near 
Seattle, Washington; traversed parts of Lake Michigan and 
Lake Erie and proceeded into the Atlantic Ocean. Figures 
12, 13, 14, and 15 show the integration points of the S-194 
main beam mapped onto selected target areas in Pass 14. In 
order to display the output of MAPS with the minimum of 
clutter, the program was altered to plot points in azimuth 
every thirty degrees for the first five degrees in angle off 
nadir. Figure 13 represents the beam pattern mapped on the 
earth's surface at a nadir point 42.334 degrees North 
Latitude and 83.186 West Longitude. The Great Lakes — 
Michigan, Huron, Erie, and Ontario — are represented to 
scale. It can easily be seen that the true contour of the 
northwest shoreline of Lake Erie and the southern shoreline 


of Lake Huron must be accurately represented for the output 
of ANTTEMP to be consistent with the measured antenna 





Figure 12 

























Figure 15 















tlemperature . Figure 14 represents the same scene but with 
the beam pattern advanced along the ground track one and 
one-half beam-widths to 41.512 North Latitude and 81.426 
West Longitude. Here only the contour of Lake Erie is truly 
significant to the input for ANTTEMP . Figure 15 shows the 
scene with the beam pattern one beam-width advanced from 
Figure 14. Note that now it is the southeast shoreline of 
Lake Erie that is most critical. Figure 16 incorporates a 

i • 

map of the Chesapeake Bay region. We found it was. necessary 
to make this map quite detailed, to include, for instance, 
the outlets and major portions of the Potomac, Choptank, 
Rappahannock, and James River, in order to correctly model 
the area. Similar complex bay areas would have to be modeled 
with equal precision. 

The antenna temperatures calculated for Pass 14 using 
MAPS at the four points modeled in Figures 12, 13, 14, and 
15, as well as for 16 other points, are shown in Figure 16. 

The calculated points are represented by the open circles 
and the solid line represents the actual S-194 measurements. 

An approximate topological cross-section along Pass 14 with 
an exaggerated vertical scale is shown along the bottom of 

the figure to identify the topological feature being overflown 

| 

The half-power antenna footprint on the surface is the same 
size as the Lake Michigan cross-section along the ground 
track. Since the footprint is comparable in size with Lakes 
Michigan and Erie, and with the Chesapeake Bay, a significant 
portion of the antenna beam is filled by land for all inland 
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ANTENNA TEMPERATURE *K 



Figure 16 




water areas in Pass 14. Since the brightness temperature of 
the land is 266°K, more than 150°K above the brightness 
temperature of water, the antenna temperature is significantly 
increased if even a small portion of the beam is filled by 

I ' . 

land. For example, if there had been no land in the antenna 
beam, antenna temperatures 44.8°K and 58.6°K lower than the 
minimum actually measured would have been measured for Lakes 
Michigan and Erie, respectively. The brightness temperature 
for land of 266°K was obtained from the measured antenna 
temperature west of Lake .Michigan. where the beam was completely 
filled by land. This value for land was used in all of the 
model calculations for the entire pass. Due to the amount 
pf land in the antenna beam, the effects of salinity of 
about 10 °K between the lakes and the ocean are reduced by 
about 30% to 7°K. This change due to salinity is significant 
and well above the measurement error. Thus, even with the 
large land effect, the salinity difference between the lakes 
and ocean are detectable from the S-194 measurements. The 
antenna temperature measurements and calculations are in 
excellent agreement with a mean difference of 0.063°K and a 
standard deviation of 1.8°K. This indicates that accurate 
model calculations can be made even in complex land-water 
regions where the antenna temperature is changing rapidly . 

Since the third objective of this investigation is to 
measure the effects of salinity, it was hoped that a study 
of Pass 7 from the North Atlantic to the South Atlantic, 
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passing within two degrees of longitude of the Amazon River 
outlet, would show the change in ocean salinity due to the 
fresh water outflow from the Amazon River. Unfortunately, 

■| ■ i • i 

siin elevation angles were close to 90° during the critical 
part of Pass 7 and the measured antenna temperature is not 
reliable because of the uncertainty in the large possible 
correction to it by sun glint.. Salinity variations of from 
at least 10 to 15 PPT near the mouth of the Amazon River to 


36 PPT in the Atlantic Ocean are to be expected during the 
rainy season at this time of year. This salinity variation 
corresponds to an expected antenna temperature increase of 
as much as 20°K. Although antenna temperatures of 119 °K in 
.the Amazon mouth and 95°K in the ocean, an increase of 24 °K, 
were measured, this increase icannot definitely be regarded 
as resulting from salinity effects due to the uncertainty in 
the sun glint correction. 

Since salinity changes of up to 5 PPT do occur from 
point to point over the open ocean, it is possible to look 
for the salinity effect in the ocean data. In order to do 
this, it is necessary to remove the effects due to variations 
in sea temperature and marine wind speed along the ground 
track. The dependence of antenna temperature on salinity 
Was computed using all of the available ocean data. The 
data were normalized using the model dependence for sea 
temperature and wind speed, to an ocean temperature of 25 °C 


and a wind speed of 5 knots and are shown in Figure 17. The 
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dots are the normalized data points and the solid line 
represents the theoretical dependence of antenna temperature 
oh salinity. Although most of the data are clustered about 
3|6 PPT , and there is only a small range in salinity, the 
measurements are quite consistent with the theoretical 
model. 

The fourth objective is to examine the effects of 
changing thermodynamic water temperature on the antenna 
temperature measurements. In order to perform this analysis, 
it was first necessary to remove the influence of changing 
wind speed and sea surface salinity by normalizing all data 
points to a wind speed of 5 knots and a salinity of 36 PPT. 

All of the ocean measurements are used. Figure 18 is a plot 

*• 

of antenna temperature vs. sea temperature where the dots 
are the normalized measurements and the solid line represents 
the theoretical dependence. There is a large change in sea 
temperature., ranging from near 0°C to nearly 30°C. Only a 
small change in antenna temperature is expected even for 
this 30°C change in sea temperature. The measurements 
confirm this result and are in good agreement with the 
theoretical model. S-194 is clearly not ttie instrument ; to be 
chosen for remote determination of sea temperature. 

The fifth objective is to consider the effects of wind 
speed, sea surface roughness, and foam coverage on the 
measured antenna temperature. On September 2, 1973, during 
Pass 23, Skylab passed over a portion of Tropical Storm 





Christine. The excellent ground truth available for this 
storm makes it a good target of opportunity for the study of 
the effects of high winds and ocean roughness. Figure 19 
shows the ground track of Pass 23, through an outline 
representation of the storm area, with a graph of measured 
antenna temperature superimposed on it. The solid outline 
encompasses the center of the storm and heaviest cloud 
cover. The dotted line shows the extent of associated, but 
lighter, cloud cover. As can b- seen, the S-194 footprint 
never passes through the most i tense part of the storm. In 
addition, the effect of high wi:.ds is somewhat hidden by the 
circumstance that the S-194 footprint was moving from an 

area of 32 PPT salinity at the left end of the track to an 

• 

area of 37 PPT salinity on the other side of Christine as 
can be seen in Table 2. This change in salinity results in 
a 5 . 7 0 K decrease in the calculated antenna temperature in 
going from one side of Christine to the other. This is 
accurately reflected in the S-194 measurements. The increase 
in antenna temperature in the vicinity of Christine is about 
3 to 5°K. This increase may .result not only from roughness 
effects but also may in part be due to very heavy rain and 
local decreases in surface and near-surface salinity resulting 
from the rain. Thus, the contribution due to wind speed 
effects is difficult to assess. The increase of 3 to 5°K is 
consistent with a wind speed increase of from 20 to 30 knots 
above that outside the storm area. This is consistent with 
the position of the ground track along the outer edge of the 
storm. 
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DEGREES NORTH LATITUDE 




DEGREES WEST LONGITUDE 


Figure 19 



During Pass 79, over the North Atlantic, wind speeds of 
from 6 knots up to 48 knots were reported by FNWC. Since 
the salinity variation is only ±0.5 PPT over the ground 
track from 31.1 North Latitude, 69.4 West Longitude to 48.7 
North Latitude, 30.4 West Longitude, the effects of high 
wind speed on antenna temperature is more easily discerned 
than during Pass 23. There is an unmistakeable antenna 
temperature increase of nearly 7°K corresponding to the wind 
speed increase of 40 knots. This dependence of about 
0.17°K/knot is very close to the value of 0.16°K/knot expected 
as discussed earlier. 

In order to more fully investigate the wind speed 
effect, all of the ocean data where the wind speed was well 
known, i.e. where no asterisk appears in Table 2, and no sun 
glint was present, were adjusted to remove the effects of 
variations in salinity and sea surface temperature. Figure 
20 is a plot of antehna temperature versus wind speed for 
these data where all of the data have been normalized to 
25°C sea temperature and 36 PPT salinity using the theoretical 
dependence on sea temperature and salinity. The solid line 
in Figure 20 is the straight line which best fits the data 
as determined by linear regression analysis. The correlation 
coefficient of 0,87 indicates a very significant wind speed 
dependence. This is also apparent from the good fit of the 
data to the regression line over the considerable range in 
win<d speed of from 3 to 48 knots and corresponding change of 
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about 15°K in antenna temperature. The adjustment of the 

data to 25°C sea temperature and 36 PPT salinity, required 

only relatively small changes in antenna temperature. The 

total range in sea temperature over the data of from 3 to 

29°C required less than a 2°K adjustment, and the salinity 

range of from 32 to 37 PPT required less than a 6°K adjustment 

in the antenna temperature. Further, any possible atmospheric 

differences over the data could result in variations of at 

most a few tenths of a degree. Therefore, the preponderant 

effect expected is that due to ocean roughness, and it 

is only reasonable to accept it as the effect which is displayed 

in Figure 20. The wind speed dependence determined from the 

regression line of 0.155°K/knot is extremely, and probably 
% 

somewhat fortuitously, close to the value of 0.16°K/knot 
determined prior to this analysis from other data as discussed 
earlier. The present result validates the earlier wind 
speed dependence which was used in ANTTEMP, at least at 1.4 

GHz. 

A final treatment of the data was performed in order to 
evaluate the absolute accuracy of S-194 as well as the 
stability and repeatability of the measurements from one 
EREP Pass to another and one mission to the next. A difference 

j . • ' ' • ■ ■■■•'.' ■ ■ ■ • 

distribution of all the ocean data was computed by subtracting 

. j 

each calculated antenna temperature from the corresponding 
measured antenna temperature. The antenna temperature was 
calculated using ANTTEMP and the theoretical model for sea 


temperature, salinity, wind speed, and atmospheric effects. 

The plot of this difference distribution in °K versus Skylab 
Pass Number is shown in Figure 21. Although each individual 
pass may have a noticeable bias ranging from + 2.1°K for 
Pass 38 to -1.9°K for Pass 40, the mean for all the data 
analyzed is -0.0035°K with a standard deviation of 1.3°K. A 
" t" test was run and the mean value of this distribution was 
determined to be not significantly (t = -0.026) different 
from zero. This indicates excellent agreement between the 
model calculations and the average absolute level of the 
measurements. However, there are apparent small systematic 
changes in the absolute level of as much as 2°K from pass to 
pass. The standard deviation of 1.3°K includes the effect 
of this systematic change in level as well as errors in the 
determination of the environmental parameters,' instrumental 
variations and noise and errors due to data quantization. 

It indicates an absolute accuracy of about 1 . 3°K for the S-194 
sensor. 

CONCLUSIONS: 

In conclusion the results of the amalysis of the S-194 
data demonstrate that it is indeed possible to make accurate 
measurements of the earth's ocean areas from a satellite 
with a passive microwave radiometer. Although S-194, at a 
frequency of 1.4 GHz, is only relatively weakly dependent 
upon marine environmental conditions, except for its dependence 
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upon salinity, compared to higher frequency radiometers, the 
rras absolute accuracy, of 1.3°K, of the measurements over a 
wide range in environmental conditions is quite good. 
Unfortunately, the footprint of S-194 is too large for a 
detailed study of land-water areas. However, we have shown 
that it is possible to accurately model complex regions such 
as lakes and bays and obtain excellent agreement between 
calculated and measured antenna temperatures. S-194 is very 
insensitive to ocean temperature, and it was not possible to 
obtain conclusive data regarding the measurement of the 
temperature of fresh water, for which it is more sensitive, 
because of the large footprint compared to fresh water 
areas. Therefore, S-194 is of very limited, if any, use in 
measuring surface water temperature. The large footprint of 
S-194 again limited the measurement of salinity in lakes or 
bay areas because of the large correction for the portion of 
the beam filled by land. The analysis of the ocean data 
shows that salinity measurements away from the influence of 
land can be made to an accuracy of about 2 PPT. A very 
clear but weak ocean roughness (wind speed) dependence was 
found of about 0 . 16°K/knot . This indicates that S-194 may 
be used to determine marine wind speeds to an accuracy of 
about ±8 knots* In order to measure either salinity or wind 
speed to these accuracies, ground truth on the other environ- 
mental parameters must be available. Under certain homogeneous 
conditions, it may be possible to measure relative changes 
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in salinity or wind speed to an even greater precision. 

Since S-194 is very insensitive to atmospheric water vapor, 
clouds and even rain, and is only weakly affected by molecular 
oxygen, the atmosphere serves principally to bias the measure- 
ments and S-194 is not useful as a meteorological instrument. 
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